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Nuclear corrections are presented for neutrino and electron induced reactions in a pedagogical manner. The 
formalism is demonstrated with numerical studies and is shown to produce substantial corrections in channels 
where the pions have the same charge with the exchanged current. Two comparisons with available data show 
consistency of the model. Additional experimental results along these lines will improve the accuracy of the 
predictions and enhance the discovery potential of experiments. 



1. Formulation of the Model 

It is widely recognized that nuclear corrections 
in neutrino reactions play an important role be- 
cause they alter sometimes the neutrino-nucleon 
interaction at the 10-30% level. There are stan- 
dard questions that are frequently asked on this 
topic for which the answers are not yet confirmed 
by observations. Among them I will mention 
three. 

(i) How big is the absorption of hadrons and in 
particular of pions within a nucleus? 

(ii) How does nuclear matter affect the propa- 
gation of resonances themselves? 

(iii) When a resonance decays within a nucleus, 
what are the changes in the energy and 
angular distribution of the decay products, 
due to interactions within the medium? 

Such questions are important for the interpreta- 
tion of the new experiments and the precise de- 
termination of neutrino parameters which occur 
in neutrino oscillations and the discovery of CP 
violation in the leptonic sector 1 . 

This article tries to present a consistent and 
easy to follow summary of the calculations in the 
ANP model [5] and then mentions a few com- 
parisons with data which are now possible. First 



I will discuss a general framework within which 
nuclear effects are calculated. 

Resonances like the P33, Pn, Sn and D13 are 
produced within nuclei with the cross sections 
proportional to the density of nuclear matter. 
They either propagate within the nuclear medium 
or decay producing hadrons like pions, a proton 
or a neutron, which propagate in the medium be- 
fore they exit from the nucleus. In this article we 
assume that the resonances either decay imme- 
diately or are absorbed. The latter means that 
they interact with the medium producing excited 
nuclear states which subsequently vibrate decay- 
ing eventually to the ground state. The extra 
energy in this case is transferred to thermal or 
other unobserved energy. The absorption of pi- 
ons will be parametrized by a phenomenological 
function aabs{W), with W the invariant mass of 
the pion nucleon system j^j. 

In this article we shall consider the pions pro- 
duced from the decays of the resonances and the 
interaction of the pions with an isoscalar medium. 
The propagation of the pions is characterized by 
three eigenfunctions with characteristic eigenval- 
ues A^: 
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The eigenfunctions Qi, i = 1, . . . , 3, denote the 
pion charge multiphcities (populations) after suc- 
cessive scatterings. More precisely the pion ad- 
mixtures included in the eigenfunctions are repro- 
ducing themselves after each scattering. Assum- 
ing charge-symmetry for the interactions there 
are three independent transport functions. In 
addition, assuming A-resonance dominance the 
elementary interactions are written in terms of 
the TT~^p cross section. The probability of a pion 
origninally at point y to propagate to point x and 
interact there is governed by an exponential law 
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with K = p{x)(Tg{W) and crg{W) an effective 
cross section given by 
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In this form, the effective cross section depends 
on the absorption which is assumed to be the 
same for the three types of pions plus an interac- 
tion cross section which includes charge-exchange 
terms. The transport problem is a matrix prob- 
lem which is solved as an eigenvalue problem. 
The absorption term, being proportional to the 
unit matrix, brings an overall factor 
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consisting of a Pauli suppression factor g{Q^,W) 
at the production point times a transport func- 
tion /(A = 1, W) for the eigenvalue A = 1. The 
eigenvalues and eigenfunctions in eq. (1) corre- 
spond to the interaction matrix, i.e. the second 
matrix in eq. As we describe later on, each 
eigenvalue has its own transport factor /(A, W). 



The absorption cross section is a phenomeno- 
logical function which contains several nuclear ef- 
fects; for example excitation of the nucleus, a 
short propagation of the A-resonance before it 
decays, etc. It is determined phenomenologically 
from data. 

Let us denote by dcr* with i = 0, — the 
neutrino-nucleon cross section averaged over pro- 
tons and neutrons at production point within the 
nucleus. Similarly denote by da[T^ , tt*) the num- 
ber of pions with charges i = 0, — emerging 
from the nucleus T^. The isospin symmetry of 
the problem gives the following relation 



da(T^; 7r+) 
d<j{T^- ttO) 

d<7{T^- TT~) 



= A{Q^, W) M 



da {71+) 
da{TT°) 
da{Tr^) 
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with 
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The parameters c and d depend on the charge 
exchange and are given in ref. [2]. They will be 
defined also below. Taking the sum of all pions 
we arrive at the relation 



da{T^; n') = A(Q^ W) V MiKdain'') . (7) 



When the Pauli factor is weighted by it 
gives a value between 0.96 and 0.98 and conse- 
quently ^(H/^) gives practically the suppression 
introduced by absorption. This is seen in table 
n where the Pauli factor is given as a function of 
Q2 in GeV2. 





9iQ^ Ma) 


0.05 


0.87 


0.20 


0.98 


0.40 


1.00 



Table 1: Pauli factor in dependence of Q^. 
We describe next the transport functions. 

In general the transport problem is a random- 
walk problem of the pions within the nucleus. 
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which is usually included in Monte-Carlo pro- 
grams as a multi-scattering process. In ref. 
there are analytic solutions for the general case as 
well as approximate solutions for specific cases. 

Since the pion-nucleon cross section in the A- 
region is peaked in the forward-backward direc- 
tion {(7 (6) ~ 1 -I- 3cos^6') we shall consider an 
average of the cross section in the forward and 
backward hemispheres. We define a second Pauli 
factor h±(W) corresponding to the scattering of 
the pion on a bound nuclcon. The Pauli fac- 
tors h±{W,cos6) are given in Appendix C, part 
2 (page 2142) of ref. [2]. At the same time we 
introduce an average over forward and backward 
hemispheres 

{a^+p{W)h+{W)) = j\^+p{W, cos d) 

X h+(W, cos 0)d cos (8) 

with a similar integration defining the averaging 
over h_{W). Then the probability for transmis- 
sion in the forward direction is given by 
1 rL{b) 

U{X^L,W)^— d2/e-^^(i"^^+) (9) 



with 

yK{l - A/X+) 



yp{Q)[aabs{W) 
1 



\{a,^^XW)h+{W))] ,(10) 



A one of the three eigenvalues, and L{h) the ef- 
fective profile of the nucleus at a given impact 
parameter h. In other words, when the neutrino 
produces a pion at an impact parameter 6, the 
pions begin to rescatter forward along the line of 
the impact parameter until they exit from the nu- 
cleus. For each line with impact parameter h the 
effective density profile is given by 



L{h) = ^/^R 



+ 




(11) 



where R is the radius of the nucleus and the den- 
sity is parametrized as 



p{r) — p{0) e 'n^ 
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given in the Landoldt-Bornstein Tables PJ- 

When the scattering is in both forward and 
backward directions the solution is known and 
given in eqs. A25-A27 of ref. 12. They define 
the transition probabilities /±(A) corresponding 
to the two directions. In terms of them we define 
_ 1 1 /(5/6) 1 /(1/2) 
3 2 /(I) 6 /(I) 



(13) 



and 
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/(I) 
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and the matrices M± in analogy to eq. (6). Then 
using differential cross sections we can obtain in- 
formation on the corrections to the original angu- 
lar distribution. To be specific we define differen- 
tial cross sections 



da{T^,q) and dcr{q) 



(15) 



for the pions emerging from the nucleus and pions 
from the production point, respectively. The unit 
vector q denotes their angular direction. Then 



da{Tq) = M+da{q) + M_dq{-q) 



(16) 



smce 



M++M-=M and /+(A)+/- (A) = /(A) .(17) 
The above equations imply 

d<7{T, q) + d(j{T, -q) = M [da{q) + d(7{-q)] .(18) 

A difference between forward and backward scat- 
tering in eq. (8) comes from the /i±(VF, cos 0)) 
functions. In this article we present results for 
the charge-exchange corrections integrated over 
angles. It will be interesting to study in the fu- 
ture the modification on the angular distribution 
caused by the nuclear rescattering j^. 

This completes the formalism for calculating 
the transport matrix in terms of the absorption 
and pion-nucleon cross sections; all other param- 
eters are defined by properties of the nuclei. 

2. Experimental Consequences and Com- 
parisons 

There are two ways for using this model. The 
first one uses experimental data in order to deter- 
mine the parameters. The second method is to 
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use nuclear parameters from tables and calculate 
the charge-exchange matrix. Since the available 
data is very limited, one uses the second method 
to compute the exchange matrix and then com- 
pare them with the few available experimental 
numbers. 




Figure 1. Electroproduction for tt" on Oxygen tar- 
get. 




Figure 2. Electroproduction for 7r+ on Oxygen tar- 
get. 



Before I present numerical values, I will state a 
general result for medium heavy nuclei. In a lep- 
ton nucleus interaction the pions which have the 
same charge as the exchanged current are reduced 
substantially, by approximately 30-40%, and for 



pions with different charge there is a slight in- 
crease. For instance, in electroproduction and 
neutral current reactions the large reduction oc- 
curs for TT^'s, while for charged-current neutrino- 
induced reactions the large correction occurs for 

TT+'S. 

Figures n and 13 show the electroproduction of 
pions on sO^^- The three curves were computed 
as follows: the dotted line is the cross section 
without nuclear corrections (average over protons 
and neutrons), the broken-line curve is the cross 
section within the averaging approximation [2] 
and the solid curve is the cross section including 
nuclear corrections without making any averag- 
ing, i.e., the exact transport problem. We notice 
that in fig. the 7r° section is largely reduced (up 
to 40%) and in fig. El the 7r+ is increased by 5%. 
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Figure 3. Charged-current neutrino production of 
tt'^'s on Oxygen target. The solid line takes into 
account all nuclear corrections and dotted line rep- 
resent without nuclear corrections. 



Figures 13 and 01 show the same results for 
charged-current neutrino-production where the 
TT^ is largely reduced and the 7r° yield slightly 
increased. In Table El we present the transport 
function /(A = 1) for several nuclei where the ef- 
fect of the absorption cross section is now evident. 
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Figure 4. The same as in Fig. |2]for tt+'s on Oxygen 
target. 





/(A - 1) 




0.81 




0.80 


18A« 


0.65 




0.63 



Table 2: The transport function /(A = 1). 

In addition, we computed the charge-exchange 
matrix for several values of the absorption cross 
section and found that it affects A by noticeable 
amounts but the 3x3 charge-exchange matrix 
very little ( 3%). Values for the charge-exchange 
matrix were reported already and we give values 
for several nuclei. 



Carbon 



.4 = 0.79, Af(6Ci2)=A 



Oxygen 



A = 0.78, M{sO^^')^A 



( 0.83 0.14 0.04 
0.14 0.73 0.14 
\ 0.04 0.14 0.83 



/ 0.78 
0.16 
\ 0.06 



0.16 0.06 
0.68 0.16 
0.16 0.78 



Argon 



A = 0.66 , MdsAr"*") = A 



I 0.73 
0.19 
y 0.08 



0.19 0.08 
0.63 0.18 
0.18 0.73 



These matrices have been used for analyzing ra- 
tios of neutral to charged currents j^I- We men- 
tion two examples in order to show the changes 
that are introduced by rescatterings on nuclei. 

Example 1 

This case considers the ratio of neutral to 
charged-current reactions 



(19) 



which has been measured in two neutrino experi- 
ments |6I7| . The Columbia-BNL experiment had 
a complex target consisting 75% Al and 25% C. 
They reported the value [H] 



E'daAr ) = 0.17 ±0.06. 



(20) 



The Gargamelle experiment used a bubble cham- 
ber filled with Freon CFsBr and reported 



0.11 < i?'(Br) < 0.22. 



(21) 



The two experiments are consistent with each 
other. The values are, however, much smaller 
than the value calculated for the scattering on 
free protons and neutrons to be denoted by 



R{Nt) = 0.44 for sin^ 



w 



0.22. 



(22) 



The agreement is restored |S] when charge- 
exchange corrections are included 



-RthoorydsAl^^) — 0.22 



(23) 



This example shows the important role that nu- 
clear corrections played in establishing the isospin 
content of the neutral current. The same ra- 
tio could be, hopefully, measured in K2K, Su- 
perkamiokande (SK) and Miniboone. For the 
long-base-line experiment K2K there is an addi- 
tional bonus: the ratios should be larger in SK be- 
cause as the i-'^'s travel to Kamioka they oscillate 
to Vr^s which do not contribute to the charged- 
current cross section in the denominator. The 
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precise value of the ratio for the two detectors is 
calculable. 

Example 2 

Another realistic case occurred in the 
Gargamelle experiment where they measured the 
production of pions in an enriched Freon target. 
The average energy of the neutrino beam was 
low, in the A-resonance region. They observed 
the ratio 0: 



— = 2.3 ±0.9. 



(24) 



The cross sections for E^, = 2.0 GeV in units of 
10~^*cm^ are 

a{i^p^ H-pTr+) = 0.70 ±0.10, (25) 
cr(z/n ^ ^"Ptt") = 0.20 ± 0.05, (26) 
criiyn^ fi-mT+) ^ 0.20 ± 0.07. (27) 

It follows from these numbers that the ratio, av- 
eraged over protons and neutrons, is 



4.50 , 



(28) 



free - 



which is far away from the value reported by 
Gargamelle. After we include nuclear corrections 
for Bromine, the ratio becomes 



2.64 



(29) 



nuclear corrected 



in good agreement with the experimental number. 

Further information on the pion absorption 
cross section could be extracted ^ from an Ar- 
gonne experiment on Deuterium |10| and data 
from the BEBC experiment [TI] on a Neon tar- 
get. Deuterium was considered to be an isoscalar 
target of free protons and neutrons, and nuclear 
corrections were introduced for Neon. Unfortu- 
nately, the two experiments were at different en- 
ergies. However, it was possible to re-weight the 
data to the same atmospheric neutrino energy 
spectrum ^21 allowing for a direct comparison. 
The parameters of the ANP model were found to 
be in good qualitative agreement with the data 



The comparisons show that the qualitative 
properties of the pion transport model are con- 
firmed by experiments. The remaining question 
is to test the degree of its accuracy. This can 
and should be done in the running and upcoming 
experiments because the model provides a means 
for studying the mixing and other properties of 
neutrinos to a higher degree of accuracy. Precise 
determinations of nuclear corrections are also 
important for the discovery of CP-violation in 
the leptonic sector. 
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